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Abstract 

<N ■ 

<^ • The world-line (Fock-Feynman-Schwinger) representation is used 

for quarks in arbitrary (vacuum and valence gluon) field to construct 
the relativistic Hamiltonian. After averaging the Green's function 
of the white qq system over gluon fields one obtains the relativistic 
Hamiltonian, which is matrix in spin indices and contains both positive 
Qh! and negative quark energies. 

pi 1 The role of the latter is studied in the example of the heavy-light 

meson and the standard einbein technic is extended to the case of the 
matrix Hamiltonian. Comparison with the Dirac equation shows a 
good agreement of the results. For arbitrary qq system the nondiago- 
^ ■ nal matrix Hamiltonian components are calculated through hyperfine 

interaction terms. A general discussion of the role of negative energy 
components is given in conclusion. 
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1 Introduction 

The quest for the Hamiltonian which contains main features of QCD - con- 
finement and Chiral Symmetry Breaking (CSB) exists from the very be- 
ginning, when fundamental field-theoretical QCD Hamiltonians, have been 
constructed in different gauges pQ. Unfortunately (nonlocal) confinement 
cannot be seen in these local FTh Hamiltonians and for practical purposes 
another sort of Hamiltonians - Effective Hamiltonians (EH) have been mod- 
elled containing minimal relativity and string- type potentials [2]. A lot of 
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information was obtained from these Hamiltonians and the general agree- 
ment of calculated meson and baryon masses with experiment is impressive 
[3], with some exceptions for mesons (e.g. pions, scalar nonets etc.) and for 
baryons (the Roper resonance and his companions, A(1405) etc.). 

The two main defects of effective Hamiltonians are that: i) The clear-cut 
derivation from the basic QCD Lagrangian was absent and therefore it is not 
clear what are approximations and how to improve EH systematically. 

ii) Connected to that, the EH contains a large number of parameters 
in addition to the minimal QCD number: current quark masses and string 
tension (or Aqcd)- Typically this additional number of parameters is more 
than ten for detailed spectrum calculation. The most important for hadron 
masses are constituent quark masses, and the overall negative constant 
Co of the order of several hundred MeV. 

With this number of arbitrary parameters the QCD dynamics in hadrons 
cannot be fully understood and one needs another approach. This approach 
(it will be called the QCD string approach (SA)) was suggested more than 
a decade ago jlj, where the SA Hamiltonian for spinless quarks was derived 
and Regge trajectories have been obtained both for mesons jl] baryons [5], 
and for glueballs Later on, the formalism was put on more rigorous 
basis in [7j and the einbein technic jH] was used in [7] to take into account 
the string moment of inertia and to obtain the correct Regge slope. In 0]- 
[7j the constituent mass was defined using the einbein technic through the 
string tension and current quark mass; the subsequent calculation of baryon 
magnetic moments [§] has confirmed the validity of this approach. 

Another mysterious problem - of large negative constant Co - was under- 
stood recently in the framework of the same QCD string approach and Co 
was identified with the large quark self-energy term [10J. The latter is cal- 
culated through the string tension and quark current masses again without 
introduction of new parameters. 

The Spin-Dependent (SD) part of the SA Hamiltonian was calculated 
earlier [TI]. It was shown that even for light quarks one can calculate the 
leading SD terms without recurring to the 1/M expansion, but using instead 
the lowest (quadratic) field-correlator approximation ^2] which works with 
accuracy of few percents [T3] • 

The final form of the SA Hamiltonian was used to calculate the masses 
of light mesons [TU Uj)], heavy quarkonia [TBI El EE], heavy-light mesons 
UJ 121 I2JJ , baryons [22 123 121 , glueballs 123 121 , hybrids [23 [23 [23 EDI , 
gluelumpsjST]. For a review see also [23 E3 EM!- It is remarkable that in 
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most cases the agreement with known experimental data and lattice data 
is good, however only the minimal QCD set of parameters was used with 
addition of standard a s {Q 2 ). 

The most important exception in mesons from the agreement above was 
for pseudoscalars (71", K, 7], rf) which need the chiral dynamics absent in the 
SA Hamiltonian. To overcome this discrepancy it was realized in "3*o"] E""f"J 
13*7] that the chiral dynamics brings a new tadpole term, which should be 
accounted for in computation of the Nambu-Goldstone spectrum. As a result 
the Gell-Mann-Oakes-Renner relation was found in 131"] and the quark 
condensate was computed [37] in terms of the SA Hamiltonian spectrum. 

This connection allows to calculate the spectrum of Nambu-Goldstone 
mesons and their radial excitations in terms of the SA Hamiltonian spectrum, 
without introducing new parameters. 

So far so good, but to proceed further one should look carefully into the 
approximations done and understand how to improve the SA Hamiltonian 
systematically. 

The systematic procedure of the derivation of the SA Hamiltonian is given 
in [I]-[7J and discussed later in [HE]-|in]- It contains three typical approxi- 
mations: 1) Replacement of the Wilson loop by the minimal area expression 
and neglect of gluon excitation, which amounts to the neglect of mixing of 
a given hadron with all its hybrid excitations. As it was shown in [3T] the 
effect of mixing is indeed small except for the cases of states almost degen- 
erate in mass. 2) The use of the local Hamiltonian which appears in the 
limit of small gluon correlation length A (denoted as T g in most previous 
publications). Since A = 0.2 fm and much smaller than typical hadron sizes, 
this limit is legitimate. 3) The use of only positive solution for the stationary 
point equations in the einbein variable, corresponding to the positive con- 
stituent quark mass. This latter approximation means neglect of the quark 
negative energy states, and it is the main point of the present investigation. 
As a result we shall obtain the Hamiltonian containing both positive and 
negative quark components, and estimate quantitatively the importance of 
the latter. 

The paper is organized as follows. In the section 2 we derive the Green 
function for the qq system and consider in section 3 the one-body self-energy 
corrections for the quark and antiquark. Having fixed that, we turn in section 
4 to the heavy-light qq interaction and derive the corresponding qq Hamil- 
tonian in the full relativistic form, containing Negative Energy Components 
(NEC) and compare numerical results of matrix Hamiltonian with those for 
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the Dirac equation. In the section 5 the effects of NEC are derived for the 
general qq system. The last section is devoted to conclusions and outlook. 
Appendix 1 is devoted to the derivation of path integral form of the FFS 
type, in particular a novel first-order form of FFS is obtained for one par- 
ticle in external nonabelian field using Weyl representation for 7 matrices. 
Appendix 2 contains details of self-energy correction. 



2 The quark- ant iquark Green function 

We recapitulate here the steps done in derivation of the SA Hamiltonian 

m en i32i 

One starts with the Fock-Feynman-Schwinger Representation (FFSR) for 
the quark (or valence gluon) Green's function in the Euclidean external glu- 
onic fields [3U , which is exact and does not contain any approximation: 



roo px 

S(x,y) = (m+Dy 1 = (m-D) ds(Dz) xy e~ K P A exp(ig A^dz^P^x, y, s) 

JO Jy 

(1) 

where K is the kinetic energy term, 

and m is the pole mass of quark, and z M (r)— the quark trajectory with end 
points x and y integrated over in (Dz) xy . 

The factor P a (x, y, s) in ((TJ) is generated by the quark spin (color-magnetic 
moment) and is equal 

Pa(x,y, s) = P F exp[g f a^F^z^dr], (3) 

where = ^(7/x7 v lulfj,): and Pp(yPj\) in (JHJ), (0) are respectively ordering 
operators of matrices F llu (A fl ) along the path z^{r). In what follows the role 
of the operator P a (x,y,s) will be crucial, and it is convenient to rewrite 
CfiuF^u in 2 x 2 notations 

fcrBvE\ 



where er are usual Pauli matrices. 
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The next step is the FFSR for the hadron Green's function, which for the 
case of the qq meson is 

/•oo />0O ^ _ A 

G q q(x, y; A) = / ds ds'(Dz) xy (Dz') xy e- K - K tr(T(m-D)W aa (x,y)T(m'-D>)) 
Jo Jo 

(5) 

where T and T = T + are 1, 7^, 75, (7^75), "tr" means trace operation both 
in Dirac and color indices, and 

W aa >(x,y) = P A exp(ig A ll dz ll )P a {x,y,s)P' IJ {x,y,s'). (6) 

JC(x,y) 

In © the closed contour C(x,y) is along trajectories of quark z^{t) and 
antiquark z' v {t'), and the ordering P4 and Pp in P ai P' a1 is universal, i.e. 
W a(T i(x,y) is the Wegner- Wilson loop with insertion of operators (@J) along 
the contour C(x,y) at the proper places. 

The FFSR Eq. (jSJ) is exact and is a functional of gluonic fields A^F^, 
which contain both perturbative and nonperturbative contributions, not spec- 
ified at this level. 

The next step, containing important approximation, is the averaging over 
gluonic fields, which yields the physical qq Green's function G qq - 

G q§ (x, y) = (G qq (x, y; A)) A . (7) 

Here the averaging is done with the usual Euclidean weight exp(— action), 
containing all gauge- fixing, ghost terms, the exact form is inessential for what 
follows. To proceed it is convenient to use the nonabelian Stokes theorem 
for the first factor on the r.h.s. in (jUJ) and to rewrite the average of (JHJ) 

as 

(W aa '(x,y)) = (expig J dTc fl „(z)F^(z)) = 
= exp £ { -^f f dTr(l)... / d-K(n)((F(l)...F(n)}) (8) 

n=l n - J J 

where 

dn^{z) = ds^ v (z) - ia^dr (9) 

and ds^ is the surface element. In (JEJ) we have used the cluster expansion 
theorem and omitted indices of dn(k), F(k) implying dii(k) = dn^^^z) = 
ds^ u (u k )-icr^ kUk dr k , F(k) = F^ kUk (u k , x ) = $(x , u k )F^ k „ k ( k )$(u k , x ), where 
y) = P exp ig P A^dz^. 
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Eq. (JSJ) is exact and therefore the r.h.s. does not depend on the choice 
of the surface, which is integrated over in ds^ u (z). To proceed one makes at 
this point the approximation, keeping only lowest (quadratic) field correlator 
((F(i)F(k))), while the surface is chosen to be the minimal area surface. As 
it was argued in [T3], using comparison with lattice data, this approximation 
(sometimes called the Gaussian Approximation (GA)) has accuracy of few 
percent. The factors (m — D) and {m! — D') in (J5J) need a special treatment 
in the process of averaging in (J7J), and as shown in Appendix 1 of [ 1 Q| one 
can use a simple replacement, 



With the insertion of the cluster expansion (JBJ) and the operator (m — D) 
from Eq. fjlOjl into the general expressions (JJJ), (JHJ), one fulfills the first step: 
the derivation of the physical qq Green's function in terms of vacuum corre- 
lators ((F(l)...F(n))). At this point it is important to discuss the separation 
of one-body (self-energy) and two-body terms in the interaction kernel (JHJ), 
together with the separation of perturbative and nonperturbative contribu- 
tions. 

3 Quark self-energy in the confining back- 
ground 

It is clear that on physical grounds it is difficult to separate out the one- 
body (self-energy) contributions for the quark connected by the string to 
the antiquark. The situation here is different in the confining QCD from 
the nonconfining QED, since in the latter electron can be isolated and its 
selfenergy is a part of the renormalized electron mass operator. For the 
bound electron in an atom, the one body and two-body contributions to the 
interaction kernel can be separated in each order in a n , as it is done e.g. in 
the Bethe-Salpeter equation. 

In case of QCD confinement cannot be excluded in any order of perturba- 
tive gluon exchanges and the separation seems to be impossible in principle. 
Nevertheless the world-line or FFS representation (JTJ), (JHJ), (JZJ) suggests a 
possible way of separating out the one-body contributions, It is based on 
the distinguishing the perimeter (L) law and area (S m i n ) law terms in the 




(10) 
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Wegner- Wilson loop, 

(W(C)) = consteM-CiL - C 2 S mm ), (11) 

where the one-body terms is associated with the coefficient C\, while the two- 
body terms - with Ci. Going over from the Wegner- Wilson loop to the qq 
Green's function and the qq Hamiltonian, the situation however is becoming 
more complicated since (G q q) is an integral over all Wegner- Wilson loops, and 
typical loops (qq trajectories) have a finite average qq separation (r) and the 
same time length T, so that both perimeter and area-law terms contribute 
terms proportional to T. At this point the FFSR is helpful, since it allows to 
separate the Lorentz-invariant self-energy (SE) terms, Am 2 , which contribute 
to the Hamiltonian (see below and in JU]) as 4^-, where 2tu = and t 
is the physical time, and s , as in (JJJ, is the proper time. One can see 
that the SE terms are multiplied by the effective length of trajectory, indeed 
AHdt = ^fdt ~ Am 2 ds. At the same time the two-body terms in the 
Hamiltonian are proportional to to (see below and in [TUI I32j). 

To calculate the SE terms explicitly we shall use the background per- 
turbation theory with the separation of nonperturbative background 
field and valence (perturbative) gluon field a M , so that the total vector 
potential can be written as 

An = + a M . (12) 

The method 03] assumes the perturbative expansion in powers of ga^, while 
B^ enters via nonperturbative field correlators known from lattice or 
analytic jSIj calculations. Accordingly we separate the contributions to the 
quark SE (we prefer to use the m 2 instead of m, since m? appears in the 
Hamiltonian both in the FFS technic and after Foldy-Wouthuyzen diagonal- 
ization of the Dirac operators) 

m 2 {fi) = m 2 pert {(i) + Am 2 p (/i) + m 2 int . (13) 

Here m 2 ert (n) is the pole mass and its connection to the MS mass is known 
to two loops (for a detailed discussion see the book 46 ) 

m£fV) = ™(™ 2 ) {l + ^a.(mP*) + 0(a 2 )} , (14) 

while the basic nonperturbative term m 2 p was found in JU] (below we shall 
find a correction to this term) and the mixed perturbative-nonperturbative 
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contribution m 2 nt is yet to be calculated. We stress that m 2 (/z) can be found 
in a gauge- invariant form only when it is computed inside the gauge-invariant 
qq or 3q Green's function, and in principle it may depend on the system where 
the quark is imbedded. 

Since we are mostly interested in the case of light quarks, the perturbative 
mass evolution is small and unimportant, and the main term appears to be 
Am^ p , which we consider now. 

Following jTU] we consider the quadratic in (crF) term in (jHJ) and expand 
the exponent to make explicit the resulting SE term (which is exponentiated 
after all, yielding additive contribution to the Hamiltonian) . One has 

(W aa ,(x,y)) S l^l + g l a ^a pX J°dr ^dr'F, u {z{r))F p x{z{r'))^ W + ..^J 

(15) 

where we have neglected the term proportional to dSp U ds Xp since it con- 
tributes to the qq potential accounted for in Hamiltonian and not one-body 
SE terms; also the mixed terms (~ dSp V (J Xp ) contribute to the spin-orbit 
potentials, also taken into account in the Hamiltonian jllj . Here Wo is the 
usual Wegner- Wilson loop without (crF) operators. The vacuum averaging 
in (JTSJ) yields in the Gaussian approximation (see Appendix of [TT] ) 

(F pu F pX W ) = {[(F^F pX ) - g 2 J ds Xf3 (F pu F af3 ) J ds yS (F pX F yS )](W )}. (16) 

Introducing scalar functions D and D\, as in |T2] (we omit for simplicity 
parallel transporters <&(u,v)) 

g 2 (Fp U {n)F pX (v)) = l{{8p P 8 vX - 5 pX 5 vp )D(u - v)+ 

+ ^[d^{h p 5p X - h x 5 up ) + d u (h x 5 pp - hpS^jD^u - v)} (17) 
with hp = Up — v p , one has 

(Tp V a pX {F^(z)F pX (z')W ) = Q[D(z - z') +D x (z- z')\- 

~ 4 / ^a/3ds a p(u)D(u - z) J cr j5 ds lX (v)D(v - z'), (18) 

where one should have in mind that in ds a p it is always implied a < (3 both 
in (|TB|) and consequently in ()18jl. 
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Using now the identities 

(Dz) xy = {Dz) xu d A u(Dz) uv d A v(Dz) vy (19) 

roc rs rr\ roc roc roc 

I ds j dn dT 2 f(s,T 1 ,T 2 ) = ds dr x / dr 2 f(s+Ti+r 2 , r l +r 2 , r 2 ) 
Jo Jo Jo Jo Jo Jo 

(20) 

one can rewrite (JHJ), (JZj) with insertion of (|T5|) as 

G q q(x,y) = G { ^{x,y)+ 

+ tr[T(m - D)A xu a fll/ d 4 uA uv a p \d 4 'vA V yr(m' - D')A yx (F fll/ (u)F p x(v)W )] 

(21) 

where we have defined 

A xu = dse- K ^(Dz) xu , K(s)=m 2 s + -j Q l-^j dr (22) 

An alternative derivation is given in [TUj . Note that (Wo) depends on tra- 
jectories entering in A zz > in (J2TJ1 . One can now take into account that when 
| a; — u\ is small, i.e. \x — u\ < T g the influence of (Wo) on A xu can be 
neglected, since (Wo) is a smooth function of its boundaries, varying when 
they are deformed at the scale larger than T g , while A xu is singular for small 
| a; — u\. Indeed in this limit neglecting the presence of (Wo) one obtains 

A (o) _ m Kjjmlx-ul) . , 

A ™ ~ 4^ \x-u\ • [Z6) 

For large \x — u\ one can use the fact, that the product of the spinless 
quark Green's function A xy and that of the spinless antiquark A yx together 
with (Wo) yield the asymptotics of the total meson mass Mo without spin 
contributions, and without self-energy corrections 

J A xy A yx (W ) ~ exp(-M |a; - y|). (24) 

Therefore we shall use for A xy at large \x — y\ the interpolation form 

. , rhKi(m\x\) ~ ~ M 

A xu {x) g u |; , rh = m + M , M « 25 
47r / |a;| 2 

(We do not need a high accuracy of (|23j) since it will enter in the small 
correction term). 
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Now inserting ([18)1 into ([21)1 one obtains the following combination. 



J( u , v) = 6[D(u -v) + D 1 (u- v)]A uv - 

— 4:<j a po 1& J ds a p{z)D{z — u) j ds 1 ${w)D{w — v). (26) 

In the first term on the r.h.s. of (J26)) A uv enters with the factors D(u — v), 
D\(u—v) which fall off with small correlation length T g . Therefore in A uv 
was taken as . Here we tend to improve this result by taking into account 
the asymptotic fall-off of A uv as in (J2~4~j) , (|23|) . This can be done replacing in 
the free propagator by the A wv from (|25|) so that the asymptotics both at 
small and large distances is reproduced. 
Identifying m 2 np from the expansion 

(m 2 +Am 2 np -D 2 )- 1 = (m 2 - D 2 )' 1 - (m 2 - D 2 )~ l Am 2 np (m 2 - D 2 )' 1 + . . . (27) 

one obtains 



Am 2 np = -J ei 4 W m ^ ( 2 ^ l) 6(D(|^|)+D 1 (|^|))+a 2 | A xu d 4 (x-u), . (28) 



We take the lattice estimate (for the quenched case [45J), Di = and the 
relation [l2j a — \ j D(z)d 2 z, and obtain 

Am 2 p = -%(t)+ °* ■ t = (m + M )T g . (29) 
vr 2(m + M ) 2 

Here <p(t) is given in Appendix 2; it is normalized as <£>(0) = 1. The first 
term on the r.h.s. of (J29|) coincides with the result ^U] when M is neglected 
(note however that coefficients before D and D\ in ^U] have a misprint, 
and should be replaced by those in (|2lj|)). while the last term in (jSHI, which 
is a correction to the first, is new. To understand the role of this term 
in creating the total mass of the meson (still without Coulomb correction 
and spin-dependent terms), we take the case of the heavy-light meson, i.e. 
when the quark is moving in the field of an infinitely heavy antiquark. The 
Hamiltonian in this case is written as in JH]-j2I] with the SE term treated 
as in PU], namely 

ffoM = m2+ 2 ^ m " P +MoM (30) 
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where Mq(u) = ^ + s(lo),s(lo) = ^J^yjs > Q (Q) = 2.338. Taking into account 
(J2HJ) the resulting expression for H q (uj) is 

W= (_^ + ^^U + M oM . (31) 
V TT 2(M + m) 2 / 2u; 

Now uj should be found from the equation jU [7\ 

dH (u) t 



duj 



(32) 



Neglecting in (|32j) the contribution of the SE term, one obtains the values 
of ujq and M , H (uj ) = M as in [4] which are shown in Table 1 all masses 
are given in GeV. 



Table 1 Mass eigenvalues (in GeV) according to Eqs. (J31|) . 
(J32|) for different values of a s 





UJQ 


e 


M 


<f(t) 




M = M + AM S£ 





0.448 


0.735 


0.96 


0.416 


-0.106 


0.854 


0.3 


0.546 


0.498 


0.771 


0.5 


-0.105 


0.666 


0.39 


0.594 


0.407 


0.704 


0.525 


-0.101 


0.603 



For a s > the values e(u), M = % +e{uj) have been calculated in |lj taking 
the color Coulomb term — into account, while AMg^ = — ^■</ , (^) an d 
t = T g M , T g = 1 GeV -1 . One can see that AM^^ is rather stable and gives 
a correction around 15% to the total mass. The correction of the second term 
on the r.h.s. of (29) to the total Am 2 is of the order of 7% for m = , so 
that the earlier calculations made without this term in and [2*2"]-[2*4"] 

would be modified by few percent. 

As the next comparison one can take the solution of Dirac equation for 
the heavy-light meson with confining and color Coulomb term present. In 
this case the SE term is absent in the first order Dirac Hamiltonian (in 
contrast to the second-order SA Hamiltonian, obtained from FFSR). The 
results of calculations, performed in [4*71 14*5] are shown in Table 2. In this 
case to compare with the SA Hamiltonian ()31|) one should take T g in (p(t) 
equal to zero, since the linear confining potential is obtained in this limit 
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(while for T g ^ there appear corrections to the linear potential calculated 
in |1H])- Hence in (J3T|) one puts y?(t = 0) = 1 and neglects as before the 
second correction term in brackets on the r.h.s. As a result one obtains for 
er=0.16 GeV 2 (all masses are given in GeV). 



Table 2 Masses and self-energy corrections according to Eq. (J31|) 
in comparison with eigenvalues of Dirac equation. 



Ot s 





0.3 


0.39 


AM SE 


-0.227 


-0.186 


-0.171 


M = M + AM SE 


0.733 


0.585 


0.533 


M D 


0.65 


0.465 


0.401 



In the last line the Dirac eigenvalues from j^7| HB| are given to be com- 
pared with the eigenvalues of (J31j) in the next line. One can see that Dirac 
eigenvalues are about 100 MeV lower. This difference can be attributed to 
the fact that in Dirac equation both positive and negative eigenvalues (the 
latter corresponding to the backward-in-time motion of quark) are taken into 
account, while in (j3Tj) only positive values of u are considered. In the next 
section we shall discuss how the negative modes (negative solutions for u ) 
can be included in the SA Hamiltonian. 



4 The matrix Hamiltonian for the heavy-light 

qq system 

We start with the Hamiltonian for the free Dirac particle H = m(3 + ap, 
which can be diagonalized using Foldy-Wouthuyzen (FW) procedure 



H = U^U, ^(j^'^-A. ,33, 

As it is explained in detail in Appendix 1 the free Green's function can 
be written in terms of as 

S(t) = ifiU ( 6 Q {t) % H) ) e~^U\ (34) 
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In the Appendix 1 also the case of the Weyl representation is discussed 
for the Dirac particle in the external fields, which gives a representation 
similar to (|34p. i.e. having the diagonal Hamiltonian of the form of Hd in the 
exponent for the time-dependent Green's function. 

We now turn to the FFS form of the quark Green's function in the external 
gluonic field, written with the help of the einbein function to |3J [7j (this 
function was previously denoted as \i in most papers) 



S q {x,y) = J Du(D A z)^ y e V J " . (35) 

After vacuum averaging this function can be associated with the Green func- 
tion of the heavy-light meson. 

Here Du is the path integration over functions uj(t), which in our formal- 
ism jU [7j is calculated by the stationary point (steepest descent) method, 
after going over to the Hamiltonian form instead of the Lagrangian path 
integral form of (|35jh namely, 



S g (x,y) = (x J Due l fo 



y), (36) 



where we have changed from the Euclidean time t to the Minkowskian time 
tu = —it, and 

^ = £ + ! + £ + < 37 > 

Solving the Schroedinger-type equation 

(la; + ar ) Vn = £n ( £l ^ n ' ( 38 ) 

one obtains 

CT 2/3 

£ »H = (2^1/3°" ( 39 ) 

where a n , n = 0, 1, 2, ... is the set of zero of Eiry functions, ao = 2.338. As a 
result for L = one has the eigenvalues E^\uj) of H (uj), equal to 
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In our previous calculations of the spectrum the stationary point in the in- 
tegration over D u was taken at the positive solution of the equation 

M| »=0, 4°>>0. (41) 
duo lw =^ n v ' 

However there is a negative solution, at uj n = —uj^ , which was neglected 
in all previous calculations. 

The Hamiltonian (J37|) refers actually to the gauge-invariant system of a 
heavy-light qQ system, where the infinitely heavy quark Q propagates along 
the time axis and is situated at the spacial origin. 

In line with the Hamiltonian (}3*3*j) . we can write the total Hamiltonian of 
the qQ system with the lower Dirac components as 




h+- 
-h 



,ho<p = E®<p (42) 



In h the spin-dependent term g ^ ' in the exponent of (J5fi|l does not con- 



tribute to the diagonal part of (|42|) for s-wave states of heavy-light mesons, 
except for the diagonal SE term considered in the previous section (this term 
can be added replacing m 2 in (|37|). (|4T7|) by m 2 + Am^). 

Now we turn to the calculation of the terms h + _ = h*_ + . To this end 
one can use Eqs.(|5|). ©, © having in mind that for the heavy-light meson 
the Green's function of the heavy quark reduces to the parallel transporter 
y) along the straight line and the factors ds\Dz') xy e~ K and absent. 
Having in mind (j3J) one must calculate the factor 



cxp 



g 2 



d7r(l)d7r(2)(F(l)F(2)) 



(43) 



In the product d7r(l)d7r(2) the term ds(l)ds(2) contributes to the linear in- 
teraction and is present in H (uj), the term dr{l)dr(2) was calculated in the 
previous section and in [TU] and was taken into account in Am^. The mixed 
terms can be written as 

exp \-g 2 f d Sflu {u)dT{F^{u) ( ^ ' ' ' ' (44) 



z(r)l 



The diagonal terms in (j44j) contribute to the spin-orbit interaction, computed 
in ^Tj, and vanish for the s-states, while the nondiagonal terms contribute 



14 



to h + -, h_ + and calculated below. Writing 



u 



ds^F^iu) = dsi±Ei(u) + ds ik F ik = n^d uE^u) + dsB, n = — , (45) 



u 



one should average this term with nondiagonal component crE and neglect 
the last term, since the correlator (|4l)|) (BiE^) is proportional to |^ and 
small, as a result one obtains the integral 

J d 2 u((nE(u))(a-E(z(T))) = <rn J d 2 uD(u - z) = ana (46) 

where we have taken into account the definition a — \ J d 2 uD(u), and the 
fact that z(t) lies on the quark trajectory, which is the boundary of the 
integration surface. 

Finally one must replace dr in f)44j) by dt, having in mind that upper 
matrix elements refer to the positive time evolution, while lower ones (corre- 
sponding to the negative energy eigenvalues) refer to the negative (backward) 
in time evolution, viz. 

*(; ;)-(-*-«)■ 

As a consequence from (|44"|) and (f47)|) one obtains 

%(7 %(T 

/ i+ _ = -( CT „), IL + = --(«ra). (48) 

The energy eigenvalues of the matrix Hamiltonian (42) are obtained in 
the usual way from the equation 

^(h~ e - h C- E )= ' *-^<«> w 



which yields 



E = ± H + (B' (50) 

where u should be found from the condition ^\ W=UJ0 = 0, which replaces the 
old condition (|4ip. and can be rewritten as 

■2 



2hoh' - 7^3 = 0. (51) 
Zuj 
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Writing h = + M (cu) , with A = ^-(f{t), one can rewrite (JSU) for 

the case m = 0, T g = 0, <£> = 1, 

^-*{i-»(s-^)^)T (52) 

where Cl = f , w = (f) 3 ' 4 , A = ^g^. 
Solving (|52|l one obtains for cr = 0.18 GeV 2 

/io(w ) « 0.56 GeV, tu « 0.21 GeV (53) 

and the energy eigenvalue (|5U|l equal to 

£ = E(lu ) = ±0.70 GeV (54) 

which should be compared with the Dirac eigenvalue from Table 2, Ed = 
l.&lQ^fa = 0.686 GeV. Thus taking into account matrix structure of the 
Hamiltonian diminishes eigenvalue by ~ 0.1 GeV and yields values in good 
agreement with independent calculation of Dirac equation. 

It is of interest to compare this Hamiltonian with the Hamiltonian of 
the Dirac equation for a light quark in the static source of linear potential, 
considered in [J7J US] • 

H D i r ac = ap + f3(m + err). (55) 
For the solution ip n {v) represented in the form 07] 

*-">^£) <«> 

the equation H.Dirac^n — tnipn assumes the form 



+ T.G n -(e n + m + ar)F n = 



f dGn 

\ f^ + *F n -(£ n -m-crr)G n = 

where «(;', I) = (j + |)s#n(Z - j). 

Equations (JSTj) are invariant under the substitution (e n , G„, F n , «) <-> 
(— £ n , F n , G n , — k). This means that for every solution with e n > 0, and 
k having the form (|56j) there exists another solution of the form 

**«-K5SSs:) (58) 



16 



which has eigenvalue — e n , —k. 

Following the idea of the FW transformation leading to (|33|) one can also 
assume that the Hamiltonian ()55|) can be diagonalized to the form 

Hoirac — ► U + f ^ ? ] U, 

V -h{-K) J 

where h(K)(p* = e n (p* and h(-K)ip~ K = e n ip~ K . 

This brings us to the eigenvalue matrix P^j) . One special feature of this 
representation is that the states ipe n and ip- En have different parities. 



5 Negative energy states for the qq mesons 

We are now considering the qq meson states made of light quarks. The Hamil- 
tonian for positive energy states was used repeatedly (see [32J for details) and 
has the form 

= + + + 2o +ar (59) 

where £j = ^nj- 7 A i = ^^f{U)-, U = (mi + M )T g . Turning now to the spin- 
dependent term (crF) in (35) and (jSJ), one remarks that in addition to the 
one-quark corrections considered in the previous section one has also the spin- 
spin term, previously treated in JI] and yielding the hyperfine interaction, 
namely the term V4. However in the derivation only the diagonal components 
of the matrix {(y^F^ v a^F a p) have been taken into account and now we shall 
look carefully into the nondiagonal terms. 

From (|SJ) one has the following contribution (Note that for antiquark the 
spin operator a$ enters with the spin opposite to a$) 




exp 



and one should replace as usual dri = ±y^-. 




z(-ri) \ ~ / z(t2)/ 

(60) 



As a result on obtains the following spin-spin terms in the Hamiltonian: 
i) from the product of diagonal components (cr^Bcr^B) one has the usual 
hyperfine interaction ^TJ |32] 



3D(r, u) + 3D 1 (r, v) + 2r 



2 ^1 
dr 2 



(61) 
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where D(r,v) = D(\/r 2 + u 2 ), and r = zi(£) — z 2 (t) is the quark- ant iquark 
distance, and we have used (fTTj) to calculate (Bi(u)Bk(v)) . Eq. (JfiTj) contains 
both perturbative and nonperturbative contributions, the latter have been 
calculated in JT] and found to be much smaller than the perturbative ones, 
which can be easily calculated using the lowest order form of Di ^1] 

Dfrt 16^ 2 



which gives the standard result 

yidiag) = Sna^a^jr) 

The matrix element of (|63j) can be written as 

a y ifa9)) _ 2a s (<rU<TW) _ Aa s ~a ( -3, 5 = \ 

{Vnf } ~ 9^ Kn[U) ~ 9(wi + wa) V +1, 5 = 1 J lMj 

where 5" = cr + |a s (r~ 2 ), for more details see Appendix 3 of [3"T] . 

ii) We now turn to the product of nondiagonal terms in (|60|). which can 
be written similarly to (|61|) as 



r(l)o-(2) 

// r V ) — To 



y(nond) (r) = - - . ^ 



■DC 



:3C 



2 , 2x^1 



3D + 3A + (3z/ 2 + r 2 ) 



<9r 2 

Keeping again only perturbative contribution, one easily obtains 



(65) 



y(nond) = _y(diag)_ ^ 

Consider now the total Hamiltonian 

H = (ff K V2) + vtT 9) )UU + n ( r d) (75)i(75) 2 . (67) 

The energy eigenvalues can be found in the same way, as it was done in the 
previous section 



E iw) = ±M Ho (uj 1 ,uj 2 ) + Vif a9) ) 2 + (vtT nd) ) 2 ■ (68) 



hf ) T\ v hf 

To illustrate the general result (|68|) we take the case of massless quarks, 
and write the energy as 



£M = ± 1 /(M W ) + ^ 2 + | (69) 
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where h (ou) is the eigenvalue of H (u>i,uj2),^i = ^2 



h (uj) 



5 

-+UJ + 



c 



1/3 



; c = a 



2/3 a(n),a{0) = 2.338. 



(70) 



00 



Also we have defined 



2a s a ( -3, S = 



9 { 1, S=l 



) 



(71) 



The crucial point is now that u is to be found as before from the minimum 
of E(u) (we assume that minimization of the eigenvalue E(u) instead of the 
operator Hamiltonian H{ui,U2) brings about a small correction as it was in 
the case of the one-channel Hamiltonian, see jjj and numerical analysis in the 
second ref. of [Zj). This section served as an illustration of positive-negative 
state mixing due to hyperfine interaction in qq mesons. For the lack of space 
the detailed analysis of the corresponding change in the spectrum will be 
published elsewhere. 

6 Conclusions 

In this paper the systematic discussion is started of the role of negative energy 
components (NEC) for quark bound states. The NEC are automatically 
taken into account in the one-body Dirac or Bethe-Salpeter formalism. In 
the latter case however the Bethe-Salpeter wave-function contains for the 
qq meson at least eight independent components and their relative role can 
be studied only numerically. A didicated analysis of NEC was done in the 
quasipotential approximation of the Bethe-Salpeter equations jlH] and effects 
of NEC was found to be significant for the spectrum of mesons. Recently 
another approach was introduced in |qTH I%5] and developed further in |oT| lo^]. 
called the Method of Dirac Orbitals (MDO), where the quark bound state is 
expanded in a series of products of individual one-body Dirac states. In this 
case also the NEC effects are taken into account, but other approximations 
are usually done (cm. motion, higher components) which require a cross- 
check of all results and comparison to other formalisms. 

The Hamiltonian formalism and in particular the SA Hamiltonian is phys- 
ically transparent and mathematically simple, it reduces to the popular Rel- 
ativistic Quark Model (RQM) Hamiltonian in its simplest form (when string 
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motion is neglected) and therefore it is necessary to understand the role of 
NEC in the Hamiltonian form. 

This is done in the present paper using the simplest bound system - the 
heavy-light meson, where the heavy quark plays the role of external field and 
therefore results can be compared to those of Dirac equation. The comparison 
proceeds in two steps. Firstly one takes in the Hamiltonian the one-body self- 
energy terms which do not mix positive energy components and NEC. Here 
a new correction term was obtained in section 3, Eq. (41) in addition to 
the old one [TU] which gives around 7% of the total. Secondly, the NEC 
mixing appears due to the nondiagonal Hamiltonian matrix elements. The 
stationary point condition for the einbein variable u should now be applied 
to the eigenvalues of the total matrix Hamiltonian E(u) 

*=(t -t )< **#-*(«»=<>. ^ = 0. (72, 

The resulting stationary values ujq and E(uq) are in good agreement with 
the eigenvalue of the Dirac operator Ed- 

E(lu ) = E D - (73) 

This procedure justifies the stationary analysis with respect to u, since 
for h (uj) the stationary point does not exist for light quarks if the self- 
energy term — 4p^- is included in ho{u), while for E{uo) the stationary point 
uo = uj always exists. In chapter 5 the first step is done for arbitrary qq 
systems, and the matrix Hamiltonian is written down explicitly. The main 
lesson here is that NEC are mixed up by the hyperfine interaction and the 
doubly nondiagonal (for both quark and antiquark) terms can be calculated 
explicitly. The further analysis and numerous applications are relegated for 
the lack of space to future publications. However already at this stage it is 
clear that NEC are very important for the structure of the wave-functions 
and eigenvalues of mesons. This is probably even more so for baryons, where 
NEC are responsible for the correct relativistic structure of baryon wave- 
functions, which is clearly seen in the values of the 9a/ 9v ratio [53J. 

The author is grateful to A.B.Kaidalov, I.M.Narodetsky, J.A.Tjon, M.A.Trusov 
and J.Weda and members of ITEP seminar theory group for useful discus- 
sions. 

The work is supported by the Federal Program of the Russian Ministry 
of industry, Science and Technology No. 40. 052. 1.1. 1112. 
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Appendix 1 



In this appendix we shall derive several representations for the quark 
Green's function in the external nonabelian field. We start with the standard 
Fock-Feynman-Schwinger Representation (FFSR) as a warm up. To this end 
one writes first the proper-time representation in the Euclidean space-time 

S = (m + by 1 = (m- D)(m 2 - b 2 )~ l = {m-D) r dse- s ^ m2 - b2) . (A.l) 

Jo 

Now one transforms (|A.1J) to the path-integral as follows 

U\ I 'dse sD l\y) = (xle^Wl^x)^^^^-!)))^)...^^^^ 1 )^). 



(A.2) 

In flA.2|) the integration over all d 4 xi...d 4 XN-i is implied and the relation 
s = eN is used. Consider now one piece of the path in (jA.2j) and write 



ee (x n \e^-^ 2 \x n ^) = (xM^e^-^-^^(p\x n ^} = 

= e *P(xn-x„_i)-e(p M -»A M (2li|S=l))^ ^ 

Integration over d A p in ()A.3|) gives 



i„,n-l = T-T^ e — +igAx ^, Ax = X n - X n -. X . (A.4) 

Insertion of ()A.4|) in (|A.2|) finally yields 

where we have used the relation D 2 = D 2 + gcr^F^. 

Note that in FFSR (|A.5J) the exponent contains 7-martices only in the 
spin term a^F^, and moreover it commutes with 75. Therefore in the chiral 
limit (m — > 0) S is odd in 7^ irrespectively of any vacuum averaging of terms 
containing and F^ v . Hence in this form one cannot describe the effect of 
the chiral symmetry breaking and one should look for other representations 
which will be the topic of other publications. 
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We start with the case of the free quark and write the Green's function 
in the energy and the time-dependent representations S(E) and S(t) (in the 
Minkowskian space-time) 

S (E) = — ^ = -J— = — -; S(t) = r S(E)e- lEt —. 

m + tp H-E m(3 + ap-E w J-oo 2tt 

(A.6) 

Consider now the Foldy-Wouthuyzen (FW) transformation of the free Hamil- 
tonian 



where 



and 



/ cos9 — sin9 \ __, / cos9 sm9 \ „, 

U = ■ a a ' ^ • a a ' (A ' 8) 

V sin cos y / V — sin cos / 



sin 29 = ° P cos 20 = , (9 = ^arctan(^). (A.9) 

VP + TTT VP + m ^ m 

Consequently one has 

= 01/^1—1/+ S (t) = [°° d f iH i +E i u+e- iEt (A.10) 
H d -E 2tt i-oo p 2 + m 2 - £ 2 V 7 

Integrating in ()A.10|) one gets finally 

S(t) = i/3U ^ ) [/+ e -VP 2 +^ 2 l*l. (A.ll) 

Another form can be given to £(-£7) using the proper-time representation 

S(E) = = ip / e-^WP+ip^ds = ipU / e- i{Hd+ipA)s dsU + . 

m + ip Jo Jo 

(A.12) 

The form (|A.12|) is interesting since it contains the matrix Hamiltonian 
in the exponent and we shall use it now to take into account external field 
A^. The form ()A.7j) is especially convenient in the nonrelativistic case when 
|p| <C m, and then also 9 <C 1, and the FW transformation is nearly diagonal. 



22 



In the opposite case, |p| 3> m, one needs to start from another representation 
of 7-matrices, namely the Weyl representation: 

(E - p<r w - mryW)* = 0; a w = ( ^ ) , 7o ™ = ( J J ) , (A.13) 

so that the Green's function in the Weyl representation is 

S W (E) = (<r w p + m 7o - E)-\ (A.14) 
Doing the FW transformation one has similarly to (|A.7|) 

V+wiiwVw = A p = ( f^ft' -V^^ ) <A15) 



where 



cos9 w ,-sm9 w \ _ it20 n 1 "icrp 



^= 7' ; " =e-*^, ^ = ^arctan^ (A.16) 
y smt^iy, cos 6V / 2 P 

and t 2 = cr 2 is the Pauli matrix in the helicity indices. 
We now use the proper-time representation for Sw 

S w = ( x \ r dse ls{(TwV+m ^ +m) \y) (A.17) 
Jo 

where = — gA^, and split the interval (x, y) into N steps as in 
u2l).A"e = s. One has 



t(W) _ I \ ie((T wT > +m^ +iV4.)\ rr \ 
1 n,n—l — \^n\^ \^n— 1/ 



d 4 p • 



^-^- 1 )^(^( n )) e ^r(«)+^W) C /+ (e w (n)). (A.18) 

(2vr) 4 

As in (1A.3J) one can write P M = p M — gA^, and integrate over d 4 p, repre- 
senting the square root terms in H^^ through the einbein function ^(^4), 

e K VP+^~ U/i n el 2fl " 2 J (A. 19) 

one can also use identity; 

c «rn = e ia(l+pl + e - to (l-^ n) . (A.20) 
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Now the integration over d 4 p yields 

li W li = j d »n (^) 2 U w (6 w y* A ^G(n)U+(0w) (A.21) 
where we have defined the diagonal matrix G(ri) with elements 
G ++ (n) = e - A -^( x2+1 ) + ^^(Ax 4 )ii^+e A -[t(x 2 + i) + #]^ ( _ A;r4) iz^ 

= e- A -[f^ 2+1 ) + ^^(A a ; 4 )i^^+e A -[t(x 2 + i) + #]i±^^ ( _ A;r4) 

(A.22) 

and x = A:r M = (x n - rc n -i) M , while p residing in 9 W is p = /ix. 
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Appendix 2 



The function <p(t),t = m/5, denned in Eq.(29) can be written as (note 
the difference in definition here and in [5|) 

<p(t) = t / z 2 dzK l (tz)e~ z (A2.1) 
Jo 

where K\ is the McDonald function, K\{x){x — > 0) ~ -, so that for t = 
one obtains 

ip(0) = 1. (A2.2) 
For t > the integration in (|A2.1|) yields two different forms; e.g. for £ < 1, 



3t 2 , 1 + 1 + 2t 2 / A r\ o \ 

^--JT^W^ i + JT^W (A2 ' 3) 

while for t > 1 one has instead, 

^) = - (^2^)5/2 arctan( v / ^T) + ^ + _^ 2)2 - (A2.4) 
For large t one has the following limiting behaviour, 

V^) = !-g + 0(i). (A2.5) 
For small t one obtains expanding the r.h.s. of ()A2.3j) 

<p(t) = 1 + t\A - 3 In + t 4 (I - y In f) + 0(t 6 ). (A2.6) 
Some numerical values are useful in applications. 

^(0.175) = 0.88, ^(1.7) = 0.234, tp(B) = 0.052 
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